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Abstract 



Let s be a point source of light inside a simple polygon P of n vertices. A path from 

s to some point t inside P is called a diffuse refl,ection path if the turning points of the 

^O path lie in the interiors of the boundary edges of P. A diffuse reflection path is said to 

be optimal if it has the minimum number of turning points amongst all diffuse reflection 

r— I paths between s and t. A diffuse reflection path between s and t inside P is said to be 

\mJ constrained if the path is simple and it intersects only the eaves of the Euclidean shortest 

\mJ path between s and t. For computing a constrained diffuse reflection path from s to i 

cc having the minimum number of turns amongst all such paths, we present an 0{n'^2'^) time 

i__i algorithm, where e > is the number of eaves. The exponential factor in the running 

time can be removed by computing a longer path that may require two more additional 

►5^ turns per eave. The constrained diffuse reflection path computed by our algorithm has at 

^> most c times the number of turns in the optimal diffuse reflection path, where 2 < c < 4. 

cn 

^' 1 Introduction 

o 

2__\ Problems of direct visibility have been studied extensively in the last few decades (see [0]). Two 

points inside a polygon P are said to be visible (directly) if the line segment joining them lies 
totally in the interior of P. The region of P visible directly from a point light source s in P is 

rS called the visibility polygon of P from s (see Figure [I]). Efficient algorithms have been designed 

cd for computing visibility polygons under various conditions [9]. In this paper, we consider a 

problem of computing the visibility of a point from a source of light by multiple reflections 

inside an n-sided simple polygon. We assume that all edges of P can reflect light like mirrors. 
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Figure 1: The region R is directly visible from s. A ray from s reaches u after one specular 
reflection. A ray from s reaches v after two diffuse reflections. 



Note that some points of P that are not directly visible or illuminated from s, can become 
visible due to multiple reflections on the edges of P (see Figure [l]). As per the law of reflection, 
the reflection of a light ray at a point is called specular if the angle of incidence is equal to the 
angle of reflection. In the other type of reflection of light called diffuse reflection, a ray incident 
at a point of an edge e is reflected in all possible interior directions except along the edge e. 
We also assume that any ray of light incident at a vertex is absorbed and not reflected. 

Multiple reflections arise naturally in the realistic rendering of three-dimensional scenes 
[H El [71 [in] . In rendering of images by ray-tracing, light sources reachable by a small number 
of reflections through an image pixel would contribute intensely at the pixel because of limited 
loss of intensity at each stage of reflection. Computing diffuse reflection paths of light arriving 
from a light source by a small number of reflections (or turning points) is therefore an important 
problem. In this paper we focus on the polynomial time computation of certain constrained 
sub-optimal paths of multiple diffuse reflections. Optimal diffuse reflections paths can always 
be computed in exponential time [ITj. Whether the computation of optimal diffuse reflection 
paths is NP-hard remains an open question. 

We state some known results on visibility with multiple reflections. In [2], Aronov et al. 
studied the region visible from a point source inside a simple n-vertex polygon when at most 
one specular (or diffuse) reflection is permitted on the bounding edges. They established a 
tight B(n^) worst-case combinatorial complexity bound for the region visible after at most 
one reflection. They also proposed an algorithm for computing such regions in 0(n^ log n) 
time for both specular as well as diffuse reflections. Aronov et al. [1] addressed the general 
problem where at most k > 2 specular reflections are used. An upper bound of 0{'n?^) and a 
worst-case lower bound of Q{{n/k)'^^) was established on the combinatorial complexity of the 
region visible due to at most a constant number k of specular reflections. They also proposed 
an algorithm running in O(n^'^logn) time, for k > 1. Davis [5], studied several variations of 
reflection problems. 



Prasad et al. [TTj showed that the upper bound on the number of edges and vertices of 
the region visible due to at most k diffuse reflections is 0(n^l^'^'^"'"^)/^^"'"^). They designed an 
(^j^^2[(fe+i)/2]+i jQg^-j time algorithm for computing such a visible region. In [T7] they conjec- 
tured that the complexity of the region visible due to at most k diffuse reflections is Q{n'^). 
Note that this region may contain blind spots or holes (see in [IE]). Aronov et al. [3] recently 
showed that the complexity of this visible region is 0{n^). Bridging the big gap between the 
0{n^) upper bound of Aronov et al. [3], and the Q{n'^) lower bound as in |[17j, is an open 
problem. 

A diffuse reflection path drp{s, t) from s to t inside P is a path inside P from s to t such that 
the turns of the path are in the interiors of the boundary edges of P. A drp{s, t) is said to be 
optimal if it has the minimum number of reflections amongst all diffuse reflection paths between 
s and t. Recently, Ghosh et al. [10] have presented three different algorithms for computing 
sub-optimal diffuse reflection paths from s to t inside P. For constructing such a path, the first 
algorithm uses a greedy method, the second algorithm uses a transformation of a minimum link 
path, and the third algorithm uses the edge-edge visibility graph of P. The first two algorithms 
are for simple polygons, and they run in 0{n + klogn) time, where k denotes the number of 
reflections in the constructed path. The third algorithm is for both polygons with or without 
holes, and it runs in O(n^) time. The number of reflections in the path produced by the third 
algorithm can be at most three times that of an optimal diffuse reflection path. 




Figure 2: A cdrp{s,t) intersects only eaves of SP{s,t). 



Let SP{s, t) denote the Euclidean shortest path from s to t inside P. Let uv be an edge of 
SP{s,t) such that if P is cut along uv, then s and t belong to the two different subpolygons of 
P (see Figure [2]). Such an edge uv is called an eave. Note that every drp(s,t) must intersect 
all the eaves of SP{s,t). A drp{s,t) is said to be constrained if (i) it is simple, (ii) it intersects 
only the eaves of SP{s, t), and (iii) it interesects each eave exactly once. Such a path is denoted 
as a cdrp{s,t). If every drp{s,t) intersets itself or intersects a non-eave edge of SP{s,t), then 
there cannot exist a cdrp{s,t) (see Figure Isl). If a cdrp{s,t) has the minimum number of turns 
amongst all cdrp{s,t), then it is denoted as mcdrp{s,t). 




Figure 3: There is no cdrp{s,t), but there exists a drp{s,t). 

A minimum hnk path between two points s and t (denoted as mlp{s,t)) is a polygonal 
path inside P having the minimum number of turns. If all turning points of a mlp{s, t) lie 
on polygonal edges, then mlp{s,t) is a drp{s,t). Ghosh [HI |9] showed how SP{s,t) can be 
transformed into a mlp{s, t) in P. Following a similar approach as shown in the sequel, SP{s, t) 
can also be transformed to a cdrp{s, t) (if it exists). So, SP{s, t), mlp{s, t) and cdrp{s, t) belong 
to the same homotopy class [HI [T3] . In addition, they all have the same number of reversals of 
turns from s to t only at eaves. 

In the next section, we present an 0{'n?2^) time algorithm for computing a mcdrp{s,t), 
where e > is the number of eaves in SP{s, t). The exponential factor in the running time can 
be removed by computing a longer path that may require two more additional turns per eave. 
If the cdrp{s,t) computed by our algorithm has k turns, then k < 2(e + opt), where opt is the 
number of turns in an optimal drp{s, t). Since e > and opt > e, k < c ■ opt, where 2 < c < 4. 

It can be seen that the path computed by our algorithm may not be useful for c > 3 as 
there already exists an algorithm for computing a drp{s, t) with c < 3 as mentioned earlier. In 
order to ascertain that c is indeed greater than 3 for the path computed by our algorithm, we 
adopt a procedure explained later to cross-check the bound in linear time. 

We conclude the paper with remarks on the diffuse reflection diameter and constrained 
specular reflection paths. 



2 Computing a constrained diffuse reflection path 

Let SP{s,t) = {s,ui,U2, ■ ■ . ,Um,t), where ui,U2, ■ ■ ■ ,Um are vertices of P. We know that 
SP{s,t) can be computed in linear time [15]. We first present an algorithm for computing a 
mcdrp{s,t) for the special case where no edge of SP{s,t) is an eave (see Figure |4]). Later we 
allow eaves in SP{s,t). 

Without loss of generality, we assume that SP{s,t) makes a right turn at every vertex of 
the path, while traversing from s to t. So, vertices of SP{s, t) belong to the counterclockwise 




Figure 4: Every cdrp{s, t) is a convex and simple path. 

boundary of P from s to t (denoted as bdcc{s, t)). Since a cdrp{s, t) does not intersect any edge 
of SP{s,t), all turning points of the cdrp{s,t) lie on the clockwise boundary of P from s to t 
(denoted as bdc{s,t)). A path from s to t inside P is said to be convex if it makes right turns 
at each turning point on the path while traversing from s to t. We have the following lemma. 

Lemma 2.1 Every constrained diffuse reflection path from s to t is a convex and simple path. 



Proof. Let cdrp{s,t) = {s,zi, . . . ,Zp,t). Since the path is simple by definition, the next 
turning point Zj+i of Zi cannot belong to bdc{s,Zi), for all i. Since cdrp{s,t) does not intersect 
any edge of SP{s,t), Zi+i must belong to bdc{zi,t). Since each turning point Zi of cdrp{s,t) 
is an interior point of an edge of bdc{s,t), cdrp{s,t) makes a right turn at Zi, for all i. Hence 
cdrp{s, t) is a simple and convex path. EQ 



Corollary 2.2 The turning points Zi, Z2, ■ ■ ■ , Zp of cdrp{s, t) = {s = Zo, Zi, . . . , ..., Zp, Zp+i = t) 
appear in the clockwise order along bdc{s,t). 

For any point p inside P, we say that the line segment pui is a left (or, right) tangent 
from p to SP{s,t) at the vertex Ui (see Figure |4]), if both Ui-i as well as Mj+i, lie to the left 
(respectively, right) of the ray emanating from p through Mj. We have the following lemma. 

Lemma 2.3 The right and left tangents to SP{s,t) from any turning point Zi of cdrp{s,t) to 
SP{s,t) lie entirely inside the simple polygon P. 



Proof. The proof follows from Lemma 2.1 due to the convexity and simplicity of cdrp{s,t). 

a 

Let CV{s, t) be the complete visible region of P bounded by SP{s, t) and Mc(s, t) such that 
the right and left tangents from every point z of CV{s^t) to SP{s,t), lie inside P (see Figure 




Figure 5: Every cdrp(s,t) lies inside the complete visible region of P. 



pi). It follows from Lemma 2.3 that a cdrp{s,t) lies totally inside CV{s,t) with turning points 
on the polygonal edges belonging to CV{s, t) as stated in the following lemma. 

Lemma 2.4 Every cdrp{s,t) lies entirely inside CV{s,t). 

Using the following property, the region CV{s, t) can be computed in linear time by travers- 
ing the shortest path trees inside P rooted at s and t [H [9] . 

Lemma 2.5 Let u and w be the parents of a vertex v G bdcc{s,t) in the shortest path trees 
inside P rooted at s and t. The vertex v belongs to CV{s^ t) if and only if both u and w belong 
to SP{s,t) (see Figure^. 

The above discussion shows that all the turning points of a cdrp{s, t) must lie on edges of 
bdc{s,t) that also belong to CV{s,t). We refer to such edges as reflecting edges. Let us first 
identify intervals on reflecting edges that can have the first turning point of a cdrp{s,t). Let 
Ml = (mil, "^12; • • •) be the intervals visible from s on reflecting edges in clockwise order along 
bdc{s,t) (see Figure [6]). We refer to these intervals as mirrors of Mi. If t is visible from any 
point Zi on a mirror of Mi, then {s,Zi,t) is a cdrp{s,t). 

Assume that t is not visible from any mirror of Mi. We identify intervals on reflecting edges 
such that every point in any such interval is visible from some point in a mirror of Mi. Note 
that a point on a reflecting edge may be visible from points in two or more mirrors of Mi. 
For each reflecting edge, the union of such intervals gives disjoint intervals on that reflecting 
edge. Let M2 = {m2i,m22, ■ • •) be the set of these intervals or mirrors in clockwise order along 
bdc{s,t). We have the following lemmas. 

Lemma 2.6 All mirrors of M2 appear after all mirrors of Mi in clockwise order along the 
boundary of CV{s^ t) . 
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Figure 6: Mirrors Mi = (mii,mi2), M2 = (^21,^.22,^23), M3 = (mai), M4 
computed in clockwise order along i?CV^(s,t). 



(m4i, 77142) are 



Proof. Let a, 6 and c be three points on bdc{s, t) in clockwise order such that a and c are 
visible from s but 6 is not. In other words, a and c belong to two mirrors of Mi but b does not 
belong to any mirror of Mi. Assume that b is visible from a. So, no part of bdc{s,t) between 
a and b can intersect bs. If bdc{s,t) between b and c intersects 6s, then b does not belong to 
CV{s,t). If SP{s,t) intersects bs, then it also intersects cs, contradicting the assumption that 
c is visible from s. So, the second turning point of a cdrp{s, t) must be on subsequent reflecting 
edges of Mi in clockwise order. 0] 

This property can be generalized for mirrors of Mi, M2, M3, . . ., as stated in the following 
lemma (see Figure |6]) . 



Lemma 2.7 The reflecting edges containing the mirrors of Mi, for the ith turning point of 
cdrp{s,t), must occur in subsequent clockwise reflecting edges containing mirrors of Mi_i for 
the {i — l)th turning point of the cdrp{s,t). 

Lemma 2.8 No point on a mirror of Mi_2 is visible from any point on a mirror of Mi, for all 



Corollary 2.9 // every turning point on a cdrp{s, t) belongs to a distinct Mi, then this cdrp{s, t) 
is a mcdrp{s,t). 

In the following lemma, we present the upper bound on the number of mirrors in Mi for 

i> 1. 



Lemma 2.10 Let ri be the number of reflecting edges containing the mirrors of Mi for i > 1. 
Then the number of mirrors in Mi is at most |Mi| + {r2 + r^ + . . . + Vi) . 



ai2 




Figure 7: After locating a\ 
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. . . and &ii,&'i25- •• on BCV{s,t), M2 = {m2i, 17122, m23, . ■ •) is 



computed from Mi = (mn, mi2, mi3, . . .). 



Proof. Consider the first mirror rriu (see Figure [?]). Let an and fen be the endpoints of mn, 
where 611 is the subsequent clockwise point of an on the boundary of CV{s,t) (denoted as 
BCV{s,t)). Let a'li be the first point in the clockwise order on a reflecting edge of M2 that 
is visible from an. Similarly, let b[i be the last point in the clockwise order on a reflecting 
edge of M2 that is visible from bu. Observe that if the right tangent from 611 to SP{s,t) is 
extended to BCV{s,t), then it meets BCV{s,t) at b'^. The portion of bdc{a'ii,b'ii) belonging 
to BCV{s,t) is called the span of mn, and is denoted as span{a[i,b[i). In the same way, the 
span of any mirror niij can be defined. Observe that b'^, b'12, &13, • ■ ■ occur in clockwise order 
along BCV{s,t). 



From the above definitions, the mirrors of M2 formed due to reflections on mn must belong 
to the portion of span{a'ii,b'ii), weakly visible from mu. A point x G BCV{s,t) is said to 
be weakly visible from a mirror if x is visible from some point of the mirror. If any additional 
mirror is formed on span{a[i,b[i) due to reflections on mirrors mi2,mi3, . . ., then there must 
exist a point w G span{a'ii,b[i) visible from some mirror of mi2,mi3, . . ., but not visible from 
mil (see Figure [?]). This is not possible due to the fact that if bdda'n, w) intersects the segment 
aiiw (other than at w), it implies that w cannot be visible from any mirror of m,i2,m,i3, . . .. 
So, if any point w' G span{a'ii,b'ii) is visible from some mirror of mi2, "n^is, . . ., w' is also visible 
from rrin. 

Let us identify mirrors of M2 formed due to reflections on m,i2. These mirrors are formed 
on the weakly visible portion of span{a'i2,b'i2) after excluding the portion span{a'ii,b'ii). If 
a'i2 G span{a'ii,b'ii), it follows that mirrors of M2 formed due the reflection on mi2 must 
belong to the non-overlapping portion bdc{b'ii,b'i2), weakly visible from mi2. If span{a'ii,b'ii) 
and span{a'i2,b'i2) are disjoint, then mirrors of M2 formed due to reflections on m,i2 belong 
to weakly visible portion of span{a'i2,b'i2)- If ■span(a'^2) ^12) contains span{a'ii,b'ii), two weakly 



visible portions of bdc 



ai2, ctiiJ, 



and bdc{b'ii, b'12) contains mirrors of M2 formed due to reflections 



on mi2. 

For identifying mirrors of M2 formed due to reflections on rrii^, remove span{a[i,b[i) and 
span{a[2,b'i2) from span{a'i^,b[^). Mirrors of M2 formed due to reflections on mis lie on the 
remaining portions of span{a[.^,b[^), weakly visible from mis. Since there is one to one cor- 
respondence between endpoints of such disjoint portions of BCV{s, t) and the endpoints of 
^^11, ''^12, • • •, the total number of mirrors in M2 is at most |Mi| + r2. 

Using a similar analysis, it can be shown that the total number of mirrors in M3 is at most 
(|Mi| + r2) + r^. Thus, the number of mirrors in Mi is at most |Mi| + (r2 + ra + . . . + Tj). 0] 

Corollary 2.11 The total number of mirrors N in Mi, M2, . . . , M^ is at most A;(|Mi| + n). 
Let us explain how mirrors of Mi,M2, . . . , M^ can be computed by traversing BCV{s,t) 



from s to t in clockwise order using the method stated in the proof of Lemma |2.10[ We know 
that each edge of BCV{s,t), partially or totally visible from s, is a mirror of Mi. Then, 
b[i,b[2, ■ ■ ■ are computed by scanning BCV{s,t) in clockwise order. The point a'^j is either the 
the first vertex y in clockwise order among all reflecting edges of M2 or the intersection point 
of the ray drawn from an through its next vertex in SP{aii,y) with the reflecting edge of 
M2 (see Figure [7]). After locating a'n, weakly visible portions of span{a[i,b'ii) are computed 
from mil. Then, depending upon position of a'l • for j > 1 with respect to the union of 
span(aii,6ii), span{a[2,b'i2), ■ ■ ■ ,span{a[j_i,b[j_^), locate the mirrors on the weakly visible 
portions of span{a[j, b'^j) after excluding spans of previous mirrors mn, mi2, . . . , mij_i. These 
partially or totally weakly visible edges are mirrors of M2. 

Similarly, mirrors of M3 can be computed from mirrors m2i,m22, ■ ■ ■■ Repeating this process, 
mirrors of M4,M5, . . . , M^ can be computed. The process involves scanning BCV{s,t) twice 
in clockwise order, once for locating all b' stage by stage, and the other for computing weakly 
visible edges from corresponding mirrors in disjoint portions of BCV{s,t). Since two or more 
mirrors can belong to the same reflecting edge of M2, . . . , M^ and the total number of mirrors 



can be Oln^) due to Corollary 2.11, the entire process takes O(n^) time. 



Let us explain how weakly visible edges from a mirror is computed. If the weak visibility 
region of a reflecting edge containing one or more mirrors of Mj is known, then all mirrors of 
Mj+i formed on its boundary can be identified by scanning it once. Since the weak visibility 
regions of every edge of CV{s,t) can be computed using the algorithm of Hershberger [12] for 
computing visibility graphs of simple polygons, the overall time complexity of the algorithm is 

Let us state how mcdrp(t, s) = {t, z^, . . . , Z2, Zi, s) can be computed from M^, Mk_i, . . . , Mi 
satisfying Corollary 2.9 (see Figure [6]). Consider the computation of Zk- Scan BCV{s,t) from 



t to s and compute the shortest path from t to each vertex v, till a point Zk on a mirror vw of 
Mk is found to be visible from t, where Zk is the intersection point of vw and the ray drawn 
from s through the next vertex of SP{t, v). Note that SP{t, v) makes only right turns and Zk is 
a point directly visible from t. Starting from z^, a similar procedure can be adopted to locate 
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a point Zk-i on a mirror of Mk-i, directly visible from Zk- Repeating this process, all turning 
points of the mcdrp{s, t) can be computed in linear time. 

Lemma 2.12 If SP{s,t) does not have an eave, then a mcdrp{s,t) can be computed in 0{n'^) 
time. 




Figure 8: There may be two turning points in a mcdrp{s,t) corresponding to every turning 
point of an optimal drp{s,t). 



Lemma 2.13 If SP{s,t) does not have an eave, then the number of turns in a mcdrp{s,t) is 
at most twice that of an optimal drp{s,t). 



Proof. Let di e bdc{s,t) be a turning point of an optimal drp{s,t) (see Figure |8]). Traverse 
this optimal path from di to t till a turning point dj G bdc{di, t) is reached. If didj is a segment 
and it does not intersect the mcdrp{s,t), then there is no turning point of this mcdrp{s,t) on 
bdc{s,t) between di and dj. If didj is a segment and it intersects the mcdrp{s,t), then there 



cannot be three or more turning points of mcdrp{s, t) on bdc{di, dj) due to Lemma 2.8 , because 



in such a case the first and the third turning points would become visible. If didj is not a 
segment, then no turning point of the optimal drp{s, t) between di and dj belongs to bdc{di, dj). 
Then, this path between di and dj still intersects the mcdrp{s, t) exactly twice, and there can be 
at most two turning points of the mcdrp{s, t) on bdc{di, dj). In the worst case, all turning points 
of the optimal drp{s,t) lie on bdc{s,t), and every link of this path intersects the mcdrp{s,t) 
twice, and therefore, the number of turns in the mcdrp{s, t) is at most twice that of the optimal 
drp{s,t). 0] 



Assume that SP{s,t) contains one eave . Let SP{s,t) = {s,ui,U2, 



^t) 



Assume tnat bi"[s,i) contams one eave . ijct bi"[s,i) = [s,ui,U2, . . . ,Uq-i,Uq, . . . ,u.rmt) 
where Ug^iUg is the eave (see Figure |9]). Without loss of generality, we assume that SP{s,t) 
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makes a right turn (or, a left turn) at every vertex of SP{s, Uq) (respectively, SP{uq-i, t)), while 
traversing from s to t. So, the vertices of SP{s,Ug^i) belong to bdcc{s,t), and the vertices of 



SP{uq, t) belong to hdc{s, t). Let cdrp{s, t) 
Uq-iUq. We have the following lemma. 



(s,^i, 



^r—l^n 



•^pi 



t), where Zj—iz^ intersects 




Figure 9: Every cdrp{s, t) lies inside the extended complete visible region of P. 



Lemma 2.14 Every cdrp{s, t) makes a right turn at zi, Z2, ■ ■ ■ , Zr-i on bdc{s, t), and left turns 



, Zp on bdc 



{s,t). 



Proof. The proof follows along the lines of the proof of Lemma 2.1 



a 



Extend Uq^iUq from Uq-i to hdc 
extend Uq-iUq from Uq to hdcduq 



{s,Uq)^ meeting it at a point Wq-i (see Figure [9|. Similarly, 
-i,t), meeting it at a point Wq. Since Zr-iZr intersects Uq-iUq, 



(i) Zr-i and Zr must be weakly visible from Uq-iUq, (ii) z^-i belongs to bdc{wq-i,Uq), and (iii) 
Zr belongs to bdcc{uq_i,Wq). Let Rq-i be the region of P bounded by bdc 



S,Wg_l), SP{S,Uq_l] 



(wq^ijUq), and Wq^iUq. Let R'q be 
Let Rq be the region of P bounded 



and Uq^iWq^i. Let R'q_i be the region of P bounded by bd^ 
the region of P bounded by bdcc{uq-i,Wq), and Uq-iWq. Ltu ±i,q 
by SP{uq,t), UqWq, and bdcdwq.t). Let ECV{s,t) be the extended complete visible region of 
U R'gU Rq such that 



Rq-l U R'q_i 



i) the left and right tangents from every point of Rq^i to 
SP{s,Uq_i) lie inside P, (ii) the left tangent from every point z of R'q_i to SP{s,Uq^i) lies 
inside P, and z is visible from some point of Uq_iUq, (iii) the left tangent from every point z of 
R' to SP{uq,t) lies inside P, and z is visible from some point of Uq^iUq, and (iv) the left and 



right tangents from every point of Rq to SP{uq,t) lie inside P. It follows from Lemma 2.14 



that every cdrp{s,t) lies entirely inside ECV{s,t) with turning points on the polygonal edges 
of ECV{s, t) as stated in the following lemma. 



Lemma 2.15 Every cdrp{s,t) lies inside ECV{s,t). 
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Figure 10: Starting from a mirror nir-i^i, 
computed across the eave Ug^iUg. 



mirrors M' 



m;i, m; 



[m 



"i,m"2) and Qi are 



As before, Mi, M2, . . . Mr-i are computed, where Mi, M2, . . . Mr-2 totally belong to reflect- 
ing edges of Rq-i, and at least one mirror mr-i,i of M^-i is partially or totally on a reflecting 
edge of R'q_i (see Figure 10). Since some points of R'g may become visible from mr-i,i, the next 



set of mirrors M^ consists of two subsets of mirrors M^ and M^', where all mirrors of M^ belong 
to reflecting edges of R'g_i, and all mirrors of M" belong to reflecting edges of R'^. Treating 
M" as Mr, the process of computing sets of mirrors Mr+i, Mr+2, ■ ■ ■ M^ continues as before till 
t becomes visible from some mirror of M^. If the process terminates before arriving at M^, i.e., 
t is not reachable, the process is again repeated treating M^ as Mr. If both processes terminate 
before arriving at M^, then there is no cdrp{s,t) in P. We have the following lemma. 

Lemma 2.16 If SP{s,t) has one eave, then a mcdrp{s,t) can be computed in 0{n^) time, and 
the number of turns in mcdrp{s, t) is at most 2 ■ opt. 

If SP{s, t) has more than one eave, the above lemma can be generalized keeping in mind 
that the process of computing sets of mirrors may be repeated 2*^ times for e eaves. 

Theorem 2.17 // SP{s,t) has e eaves for e > 0, then a mcdrp{s,t) can be computed in 
0(n^2'') time. 

The exponential factor in the running time can be removed by computing a longer cdrp{s, t) 
that may require two more additional turns per eave as follows. Let us go back to the case 



of SP{s,t) containing one eave (see Figure 10). There are two subsets of mirrors M^ e R'q-i 
and M" e R'g of Mr, and treating each of them separately as Mr, the process of computation 
is carried out twice for reaching t. These two process of computations of Mr+i from Mr can 
be combined by treating M^' U Qi as Mr, where Qi is the set of mirrors on reflecting edges 
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of R! weakly visible from mirrors of M'^. Note that since some points of M^ may be visible 
from mirrors of Qi, mirrors of M'^ must be excluded from M^+i before Mj._|_2 is computed from 
Mr+i- Thus, Ml, M2, . . . , Mj,_i, Mr, M^+i, . . . , M^ can be computed without repetition. While 
computing turning points of a cdrp{s, t) from t to s in the reverse order as stated earlier, an 
extra care should be taken to locate turning points on M^_i, Qi and M"_i. In the following, 
we present the major steps of the algorithm for computing a cdrp{s,t). 

1. Compute the Euclidean shortest path SP{s,t) from s to t. 

2. Compute the extended complete visibility region ECV{s,t) of P. 

3. Starting for s, compute mirrors of Mi, M2, . . . , M^ till t becomes visible from some mirror 
ofMfc. 

4. Starting form t, compute cdrp(t,s) = {t,Zk, ■ ■ ■ ,Z2,Zi,s) by locating turning points on 
mirrors of M^, Mfc_i, . . . , Mi. 

5. Output C(irp(s,t). 




Figure 11: An optimal drp{s,t) uses extensions of every eave of SP{s,t) as its links. For 
every turn of the optimal drp{s,t) after initial two turns, there are four turns of the cdrp{s,t). 
Observe that the cdrp{s,t) is also mcdrp{s,t). 

Let us count the number of turns k taken by the cdrp{s, t) computed by our algorithm. We 



know from Lemma 2.13 that there can be two turns in -R^-i of the cdrp{s,t) for every turn of 
an optimal drp{s,t). The same argument holds for turning points in Rg. If we assume that an 
optimal drp{s, t) uses Wg^iWg as a link in its path from s to t, then the cdrp{s, t) can have two 
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turning points in R'l and two more turning points in Rq. So, k < 2opt + 2 for one eave. If 



SP{s, t) has e eaves and an optimal drp{s, t) uses extensions of every eave (see Figure 11 ), then 
k < 2opt + 2e. Since opt > e, k < c ■ opt, where 2 < c < 4. 

Let us discuss how to cross-check the upper bound of c for the cdrp{s, t) computed by our 
algorithm. It can be seen that counting the number of additional turns actually taken by 
our computed cdrp{s, t) at each eave, a realistic tighter upper bound for c can be estimated. 
Moreover, since the number of turns in an optimal drp{s, t) is at least the number of turns of 
a mlp{s,t), the number of turns taken by our computed cdrp{s,t) can also be compared with 
that of mlp{s, t) to estimate an upper bound for c. Since a mlp{s, t) can be computed in a 
linear time [H [9] , the entire checking procedure can also be done in linear time. We summarize 
the result in the following theorem. 

Theorem 2.18 If SP{s,t) has e eaves for e > 0, then a cdrp{s,t) can be computed in Oiji^) 
time, and the number of turns in the cdrp{s,t) is at most c times opt, where 2 < c < 4. 

3 Concluding remarks 

Let a and b be two vertices of P such that the number of turns in an optimal drp{a, b) is 
maximum amongst all pairs of vertices of P. The number of turning points in such a path 
is called the diffuse reflection diameter of P. Khan et al. [H] have studied the relationship 
between the diameter of P and the number of vertices n of P. Here we establish a similar 
relationship for constrained diffuse reflection paths in P. 




Figure 12: All edges of bdc{a, b) except the first one can have turning points of a mcdrp{a, b). 
Lemma 3.1 // there exists a cdrp{a, b) in P, then the number of turns in mcdrp{a, 6) < |. 
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Proof. Let P be a spiral polygon of n vertices (see Figure 12). Let a and b be two vertices of 
P such that all vertices of bdc{a, b) are convex and all vertices of bdcda, b) are reflex. Since the 
diameter in P is between a and b, and there always exists a cdrp{a, b) in P, with all turning 
points on bdc{a, b). Observe that every edge of bdc{a, b) except the edge incident on a can have 
a turning point of mcdrp{a, b). Since the non-consecutive turning points of the mcdrp{a, b) are 



not mutually visible due to Lemma 2.8, and their visibility can only be blocked by bdcdo,^ 6), the 
number of edges of bdcda, b) must be at least one more than the number of turns of mcdrp{a, b). 
So, the bound holds for any spriral polygon P and it is tight. This scenario is the worst-case 
for any simply polygon P, where SP{a, b) has no eaves. 

Assume that SP{a, b) has an one eave Ug^iUg. Cut P into two sub-polygons using Ug^iUg. 
Since the bound for spiral polygons holds for each of these two sub-polygons, the bound also 
holds for P. If SP{s, t) has two or more eaves, then P can be cut using each eave, and the 
bound holds for the entire polygon P. ED 

Corollary 3.2 The number of turns in an optimal drp{a,b) is less than |. 

Proof. Since the number of turns in a mcdrp{a, b) is at least the number of turns in an optimal 
drp{a, b), the bound holds. ED 

For any pair of vertices s and t in P, if the angle of incidence is the same as the angle of 
reflection at every turning point of a cdrp{s,t), then the path is called a constrained specular 
reflection path from s to t (denote as csrp{s,t)). Unlike a cdrp{s,t), non-conseccutive turning 
points of a csrp{s, t) may be visible. Moreover, mirrors of {i — l)th and ith turning points of a 
csrp{s, t) may be overlapping. Following the method of depth first search for graphs, a sequence 
of reflecting edges containing turning points of a csrp{s,t) may be located if there exists one 
such path. However, this search may require exponential running time in the worst case. If the 
search is successful at an earlier stage, then a csrp{s,t) may be located in polynomial time. It 
remains open whether csrp{s, t) can always be computed in polynomial time. 
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